General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



■ASA TECniCAL NE3K»AKNlf 


»SA »-7699k 


LONG PERIODIC TTRXS IN THF SOLAR SVSTCH 
P. Bretagnon 


(IiSl-TB-7699i) LOIS PEEICDIC TFF?S II TIE M63- 1713* 

SOLiS STSTEH (latioaal iercaaatics anl Space 
Idainis^ra tioa) 26 p HC 3C3/SF 3C I CSCL 03B 

Unci as 

G3/96 Q8205 


Translation of "Tennes h longues periodes dans le svst^e solaire" 
Astronosy and Astrophysics, Vol. 30 , lo. 1 , Jan. igj*, pp. 1 * 1 . 15 *. 




V?’ --pp 


^cC£i/£^ 

'' ^cccss o£pr. 






NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
V:ASKINGT0N, D.C. 20546 JUNE 198? 


1 


OF POOR QIMLflY 


f I. J- 3 - *. 

[ MSA W_-lG®?L 

MWg'pErToDIC TERTIS IM THE SOLAR 
SYSTEM 


»- f- 


5.1u^’r982 

f r . • — - 


y. JUAc^al 

P. Bretagnon, Bureau des Longitudes 
Obsenr&toire, Baris, France 


Leo. Kanner Associates 
Redwood City, California 94063 


»«■ r - Kz:-. ' 


•S. •»-« 1, !•• 


,i:. Cl— 

[NASW-3541 

■ IS. Ty»« i* Rr^-- ;•»*►,-.»< Cr****' 


' - - ’Translation 

2 IS Asc"»«* I 

I National Aeronautics and Space Adminis- — — 

I tration, Washington, D.C. 20546 


15 . Haias 

Translation of "Termes h longues per* odes dans le syst^e 
solaire", Astroncaiy and Astrophysics, Voi JO. Ho. 1, Jan. 19?^, 
pp. lM-1^ 


(ATU-21142) 


I A At»*««cv 


The long period variations of the first eight planets 
in the solar system are studied. First, the Lagrangian 
solution is calculated and then the 

long period terms with fourth order e^entricities and inclin- 
ations are introduced into the perturbation function. A sec- 
ond approximation was made taking into account the short 
period terms' contribution, namely the perturbations of 
first order with respect to the masses. Special attention 
was paid to the determination of the integration constants. 

The relative importance of the different contributions 
is shown. It is useless for example to introduce the long 
period terms of fifth order if no account has been taken of 
the short period terms. Meanwhile, the terms that have been 
neglected would not introduce large changes in the integration 
constants. Even so, the calculation should be repeated with 
higher order short period terms and fifth order long periods. 
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SUMMARY [English language sunmary front the original text] 

He have studied the long period variations of the eight planets 
of the solar system (Pluto is excluded) . We first calculated the 
Lagrange solution. We then introduced the long period terms of 
fourth order in excentricities and inclinations in the disturbing 
function. In a second approximation we took into account the contrib- 
ution of the short period terms which provide the perturbations of 
the first order with respect to the masses. We have paid special 
attention to the problem of the determination of the integration 
cons tants. 

We began with the expansion of the disturbing function R [form- 
ula (1)]. We used the variables h = e sin n, K = cos ra, p = sin ^ 
cos fi, q = sin ^ cos fl and obtained expression (3) for the disturbing 
function and the equations of Lagrange (4). 

In the Lagrange method, one retains only the second order terms 

of the quantities h, k, p, q of the so called long period part of 

the disturbing function. The resolution of the system of differential 
equations thus obtained gives the solution of Lagrange (5) . The 
corresponding integration constants are given in Tables 2, 3, 4 and 5. 

We later introduced the long period terms of the disturbing 

function, of fourth order in the quantities h, k, p, q. These terns 

give rise to third order terms in the Eq. (6) for the variables h^, 
for example. We then substitute rumerically the Lagrange solution 

♦Numbers in the margin indicate pagination in the foreign text. 
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in these thiri order tens and hence obtain the fona (8) of the 
equation for dh^/dt. 

In a second approximation, we also irtroduced the short period 
tenss of the disturbing function. The masses are substituted nuiaer- 
ically and the terns thus found are indentical in form to those 
arising from long period terms of fourth order of the disturbing 
function and are directly added totheEq. (8). 


prgcb 

quality 


To solve the systens of Eq. (8) and (9) , we used the Xrylov- 

Bogolioubov method, which consists in seeking a solution of the 

form Cll) with a modification of the frequencies given by (12) - 

Through (12) and derivation of (11) we obtain (13). In addition, the 

substitution of (11) into (8) and (9) gives (14), so that we get tlie 

two expressions (13) and (14) for dh^/dt and dk^/dt; their third 

order parts are given in (15) by identification. It is then possible 

to determine the quantities M . N . B . and C . introduced in 

u , . , r , u , . , r , 3 3 

(11) and (12) . 


t 

I 

I 




i 
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The solutions are given by (16) and (17) and in Tables 8 to 13. 


The comparison between Tables 3 and 8 shows that the integration 

constants have been greatly modified, particularly for the planets 

Mercury, Venus, Earth and Mars. This is due to the importance of 

third order terms for these planets. Table 9 gives the modifications 

and of the frequencies as well as the new values of these 

frequencies; g. =g. +B.; s. =s. +C.. Tables 10 and 11 show the 

amplitude of the Lagrange solution calculated with the new constants; 

Tables 12 and 13 show the amplitudes M and N c of the argu- 

u,.,- u,.,e ^ 

ments of higher order. 


This work displays the relative importance of the different 
contributions: it is, for example, useless to introduce the long 

period terms of fifth order if one has not taken into account the 
short period terms. We have included the major contributions; the 
neglected terms would not introduce large modifications of the 
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constants of integration. However, the calculation should be repeated 
including long period terms of fifth order and short period terras of 
higher order. 

♦ 

Key words: planetary theory, secular ;>erturbations 

There have been several studies of long period terms in the solar 
system- Stockwell, Harzer (1895) , Hill (1897) , and more recently 
Brouwer and van Woerkcra (1950) and Anolik et al. (1969) . 

Brouwer and van Woerkom calculated the Lagrange solution for the 
eight planets and in particular investigated the Jupiter-Saturn case. 

This was a continuation of the work of Hill, who had determined a 

t 

mean perturbation function on the basis of Le Verrier's findings. /142 ' 

Brouwer and van Woerkom used this perturbation function, which had 
been extended to sixth order excentricities and inclinations, for 
Jupiter-Saturn. It is difficult, however, to determine the accuracy 
of their result because Hill empirically established some of the co- 
efficients . 

Anolik et al. dealt with the eight planet case by introducing 
all the perturbation function's long period terms up to fourth order 
excentricities and inclinations. 

Our goal was to evaluate the significance of the various long 
period terms according to their origin. We too dealt with only the 
eight planet problem. Pluto was neglected for several reasons. First 
of all, the generally accepted mass of Pluto, which previously had 
been 1/360,000 the solar mass, is now 1/1,8000,000 with a large 
uncertainty: 

= lSOOOOO± 600000. 

«'r 

Moreover, Pluto's radius vector can be less than Neptune's, with the 
result that expansions in a, the ratio of semimajor axes, of the ! 

I 

perturbation function, are no longer convergent. Finally, the intro- 
duction of Pluto's influence causes the appearance of very large 
resonances between Neptune and Pluto whose physical character is 


I 


unclear. 


ofkor«mu<v 


Lastly, we calculated the eight planet Lagrangian solution and 
then introduced the perturbation function's fourth order terns as 
well as the contribution of the hort period terms of first order 
with respect to the masses. In addition, we particulary concentrated 
on the problem of determining the integration constants because of 
the significance of the terms modifying the Lagrangian solution. 


The expansions of perturbation function R that we used are those 
constructed by Chapront at the Bureau des Longitudes. They take the 
fomr of analytical expansions in e and sin i/2, where e represents 
the excentricity and i the inclination of the orbital plane relative 
to the plane of origin. 

k= V |sin-^ j X fsin-^ j (1) 

with 

=Ji 'I -h^i . 


X being the planet's longitude, ^ the argument of the perihelion, and 

n the argument of the node. The subscript I refers to the inside 

planet and E to the outside one. The perturbation function's long 

period portion is that oart for which X and X_ are absent, i.e. in 

1 r. 


which “ ^2 ~ summation with respect to the small quantities 

e^ , Cg, sin ^E/2 done starting with zero order terms 

and then 2 , 3 , . . . 


The orbit's descriptive elements (the semimajor axis a, the 
excentricity e, the inclination i, the node argument fi, and the 
perihelion argument - are those of Newcomb. These elements are 
expressed relative to the 1850.0 ecliptic averaged over short periods, 
which will serve as our point of departure (t = 0 for 1850.0) in 
determining the integration constants of the sought after solutions. 
Also, we used more recent values for the masses of Venus, Earth, 

Mars, and Saturn than the ones Newcomb used. 

The mean motions n^ are the average observed values. The semi- 


4 


Of POOR Q0*UT» 

major axes a. are related to the values of n. by the expression 
"l^l “ constant. Now, we need for the mean motions values from 
which the secular perturbations have been removed. We therefore 
calculated the secular perturbations on the basis of Chapront's 
and Simon's work concerning the construction of planetary theory with 
secular terms. 

In the end, we used for each planet the value n of the mean 
motion defined by: 

n = - 5n 

and the value a of the semimajor axis obtained by 

n a = constant. 

We have assembled the elements adopted for the eight planets in 
Table 1. 


Table 1 

Planetary Elements for 1850.0 



a, (Vyr) 

a ("/yr) 

nj (AU) 

n- (AU) 

4* 

n 

1 

Q 

m /II 

WttfY 

S38I0I6JS93 

5 58 1025. 17.32 

0.5S709867I5 

0.587098.3460 

0.20.160.3% 

75 07 197.37 

7 0007 .00 

46 55 12724 

60(III(MN) 

Vans 

2I06MI.417I 

2l0665i.*65l 

0.7255522169 

0.7255298487 

0.00684458 

1 29 27 .34.5 

.3 25 55.26 

75 1947.41 

41 18 .HO 

Earth 1 :I95977.4496 

1 295975.6094 

1.000000021 

I.0000l81%8 

0.01677126 

lOO 21 56..30 

0 


5289(0 

Min 

689050.9354 

6S9059.2'<|7 

1.52.36914428 

1.5256791.387 

0.09526685 

5.3.3 17 52.57 

1 5102.42 

*■48 24 05.40 

509911(0 

L^nn 

109256.6.^»54 

I09;6.3.0.3053 

5.202805945 

5.202600424 

0.04825582 

1 1 .14 26.72 

1 1841.81 

98 55 58.16 

II >4 ’..3 .35 

45996.20414 

45885.621 12 

9.558845655 

9,554827.367 

0.0560(8175 

90 06 .39.55 

2 29 39.26 

112 20 51.58 

549 s 

'lianus 

l.<i426.09:S 

l5.3>-s."'5i 

19.18228185 

• 9.2171061.3 

0.0469055 

168 15 46.9 

0 46 20.54 

7.3 14 08.0 

22869 

ticpiune 

7864.698 

7845.528 

.30.05‘'542 

.30.111791 

0.0085082 

45 19 4.5.7 

1 4701.81 

1.30 08 00.2 

l“.3|4 


[Commas in tabulated material are equivalent to decimal points.] 

We chose the following variables to analyze our problem: 

/j = t'sinr7. p=sinjsinfl. (2 

l = fCosr7. g=sin 2 COsfl. 

This choice was made in order to avoid the appearance of quantifies 
expressed in e and i in the denominators of the Lagrangian eouations. 
Such quantities could cancel each other out. In addition, this is 
necessary for the resolving process because in this wav the solutions 
are expressed formally through the use of these variables and, in the 
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algorithm of solution's construction, the second members always re- 
tain the same polynomial form. 

The change in variaoles defined by (2) yields in the perturbation 
function in form (1) an expression of the form: 

I? = IS( /i'£- k) ' k'£ p'i p'itf, I/J-* cos( i , ;.| + I'l ( 3 ) 

where the summation is extended to such exponential values that 
rj^ + r 2 + + S 2 + tj^ + t 2 + Uj^ + U 2 < where u; is the order at 

which it is desired to limit the calculations. 

For the variables defined in (2) the Lagrange equations are 
written: 


Jh 

0 


■ cR 


h{l - t-'l 

1 2 

I'R 

^ 

kp 

cR 
■ ^ 

kq 

cR 

Jf'” 


na' 

ck 

r.ir 

[T7(iT: 

^i‘-T 

( /. 



2na'-i\-e')'‘ 

rq 

dk 



cR 


1(1 - 

I'R 

/»P 

cR 

(»<?_ 

PR 

dt 


na' 

ch 

»I<J'[1 -^11 

-f’l' 

-T ^ 

2h<j-( 1 - f*)‘ 

' t'P 

2na'{ I — «?')* 


^ = 


1 

c 

R 

p 

- 

PR 

- pk- - 

cR 

O- 

ph 

cR- 

dt 

T, 

!d*n - 

7 

«; 


-e'7^ 

r/. 


i'h 

2 »« 2 - 7 i -e')' *" 

ck 

--1 =- 



1 

(R 


J 

cR 

qk 

cR 

+ 

cR 

dt 


4na'{\ 


("P 

Ina'i 

1 -f-| 

' * d/. 

2nj-(l 

*' dh 

2»ifl^(l - 

^ dk 


1 _ 2 f'R 

a dt na^ c/. 


i). 2 cR i, i‘R , i'RX 1 ( rR cR\ 

Jl na ca nj'[l +(t - ' <h rk j I - e‘) ‘ \ cp ( qj 

2 2 2 
where e = h + k . 


LAGRANGIAN METHOD 


For a planet of subscript u perturbed by the seven other planets 
of subscript v, the perturbation function is written: 

4 V /i- 

r<» I >n ‘*1 

The first summation is extended to the planets inside the one under 
consideration, the second to the planets outside. We use the follow- 
ing notation; 

= a^/A (A = distance of the two planets) 


and 
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l + w. 


«>r 

I + lllf 


Limited to the second order, has the following expression; 

K = C.: t + kl + hi + k;.)-4AJi>l -t ill + pi + </f > 

+ B„{kX + ItJ'i't + SAJq,</r + p,p,) 

where C , A , B are functions of a = a /a , which is constant 
uv uv uv uv u V 

here. We thus obtain, with the notation; 


[».r] 


I + /«, 


if i>u, 


[«.«] = 


N,»H, 

I + m. 


if i<u. 


-ir = + 1 [••.1](2A^X4-B^,kx 

d( j 

% = - I Mi2Aj,, + B„hX 

Ul 

at 

-^ = + I [«.f](2^„,p.-2/l.,p,). 
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This system is written in matrix form as; 

dK 
dt 
dQ 


~=ExK, 

dt 




dt 


= -ExH, 
= -lxP, 


where H is the column vector with components h^ 


i • • ♦ r 




K 


the column vector (k„ , k„, ..., k„) , P the column vector (o„ , p„. 

Me V N * Me V 


...» Pjj) » and Q the column vector (qj^g» q^» ...» q^j) 


The sub- 
scripts Me, V, ..., N represent Mercury, Venus, ..., Neptune, 
respectively. E and I are the matrices of the linear systems in 
excentricities and inclinations respectively. 


The conventional resolution of the two Lagrangian systems gives 
the eigenvalues; 

g^, i = 1, 2, ..., 8 for the excentricities; 
s^, i = 1» 2, ..., 8 for the inclinations. 

7 
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One of the eigenvectors in the system of inclinations in zero. 
We will assume s^ = 0. 


We determine the eigenvalues A. . associated with g., and u. . 

ID 3 ID 

associated with s ^ , which gives the Lagrangian solution: 

hi= Y. 'ijMjS\n{gj( + Pj) 
j* I 
8 

^i= Z /■ijMjCos{gjt + Pj) 

J‘ > 

8 

P.= Z Pi^iV^sin(s/ + ^j) 

i 

8 

^ Z COSf.s/ + . 

J'l 


(5) 


Lastly, we calculate the 32 constants of integration , 8 ^ , 
Nj, from the values of h, k, p, q at t = 0. • 


We have assembled the eigenvalues g and s in Table 2, and 

in Table 3 we show the 32 constants of integration M, 6, N, 6. 

Lastly, Table 4 gives the amplitudes of the Lagrangian solution mul- 
8 8 

tiplied by 10 : A..M. x 10 for the excentricities , and similarly 

in Table 5, y N . x 10 for the inclinations. Frequencies g and s 
i j D 

are expressed in seconds per year. A., , y. ., M., N. are dimensionless 

ih 13 3 3 

numbers. 


Table 2 Table 3 /145 

Frequencies in " /yr Constants of Integration 


(Lagrangian 

Solution) 


(Lagrangian 

Solution) 



9 

s 

j 


/» 

,v 

Ji 


+ 5.4fii369 

- .< i'79658 

1 

0.I8I4IU4O 

87 irii:37 

0.06274851 

18 15 28.76 

• 

+ 7.3465«l 

- 6.5 "1387 

2 

0.0IW7I2 

l>»2 40 29.79 

0.00506380 

316 20 04.8) 

1 

+ 17,33 1295 

- I>.‘46:05 

3 

0.01056860 

332 56 51.67 

0.01222166 

2.54 1705.60 


+ |f!.0(M5S4 

- r.6.76lll 

4 

0.07.7(10403 

316 06 34.81 

0.02519918 

295 40 38.75 


3.7 1 1401 

0 

5 

0.04319426 

27 48 .74.76 

0.0)383974 

106 08 44,5(7 

m 

-t 22.2Sfi5.<2 

'-!I'6 

6 

004877747 

1.77 51 04.46 

(UJ0‘.>6?89 

|26 19 31.43 

1 

+ 2.70I787 

- 2.4Ck4526 

7 

(1.03140786 

106 19 25.39 

0.00S80286 

313 42 .55.41 

1 

+ 0.6.73116 

- r..^-75:o 

8 

0.0092.7780 

66 09 18.45 

0.00588386 

201 00 52.94 


(Commas in tabulated material are equivalent to decimal points.] 
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Table 4 

p 

10 . Amplitudes of Lagrangian Solution 


j 

i 

1 

2 

> 

4 

5 

6 

7 

8 ’ 

i 

Mercury 

lSl4l040 

-23ISI06 

155897 

- P0495 

2414071 

11364 

62612 

1 

72f 

Verms 

63 IS 73 

I90<)7|2 

-1275337 

1497153 

1624499 

- 55447 

61654 

nos 

Earth 

405S45 

14904t:>'' 

1056860 

- 1492703 

1624452 

24'593 

65295 

1284 

Mar> 

662S3 

264664 

3017351 

7300403 

’ 1870379 

i6ru34 

86626 

2053 

Jupiter 

703 

- 1055 

95 

55 

4319426 

- :5''3iUt) 

218379 

5968 

Saturn 

627 

1088 

750 

840 

3404356 

4^37743 

199229 

6'36' 

Uranus 

271 

265 

44 

46 

- 43S4597 

- I8I906 

3140786 

I4I05I 

Neptune 

4 

10 

3 

3 

160418 

- 13561 

- 338902 

9237x0 


y . .N . X 10®. 
ID D 


Table 5 

Amplitudes of Lagrangian Solution 



1 

2 

3 

4 

5 

6 

7 

8 

Mercury 

6274851 

- 1781583 

204668 

58171 

1383974 

13920 

- 166549 

- 72367 

Venus 

591896 

506380 

-1341594 

- 343391 

1383974 

6022 

- 95883 

- 66215 

Earth 

426404 

408232 

1222166 

226117 

1383974 

140699 

- 86614 

- 64905 

Mars 

90534 

90894 

- 1794150 

2519918 

1383974 

452917 

- 62850 

- 61461 

Jupiter 

- 1038 

655 

9 

88 

1383974 

-3l5S7g 

- 47877 

- 58494 

Saturn 

- 1328 

925 

241 

9i6 

1383974 

Te6‘99 

- 39034 

- 56ja 

Uranus 

1112 

477 

20 

86 

1383974 

- 34‘90 

8802S6 

547 IJ 

Neptune 

28 

27 

2 

9 

1383974 

- 3858 

- 103566 

588356 


INTRODUCTION OF HIGHER ORDER TERMS 

We are now going to introduce the perturbation function's long 
period terms of fourth order h, k, p, q, as well as the perturbation 
function's short order terms. 


Fourth Order Long Period Terms 


By differentiation, these terms yield third order terms, and the 
Lagrangian equation for variable h^, for example, then has the follow- 
ing form: 


Ji 




( 6 ) 


9 


ORIGINAL PAGE IS 
OF POOR QUALITY 


where is a homogeneous third degree polynomial. 

Into polynomial we substitute the Lagrangian solution (5), 
whose numerical values are given in Tables 2, 3, 4, and 5: 

• s 

i*> j«! 

• r 

Z M,J Nj sin Oj , <?, == y Sj cos 6j , 

where we have made , . = g.t + 3. and 0 s.t + 5.. Therefore only 

] 3 3 3 3 D 

the numerical values of amplitudes /« .M . and u .N. appear in this 
calculation. 


Among the values of the i and j subscripts of arguments 
(i = 1, 2, ..., 8) and (j = 1, 2, ..., 8), such a substitul^ion 
makes combinations appear in which at most only three values of sub- 
scripts i and j are involved. For example, there will be combinations 
of the type (4^^^ ^^'^5 ~ '^6^' 


The expression Z [“'OPy, therefore has the form: 

I • « 


I i-. wsi'.*', + ■■■ + i,<e, + ;,»,+ ■ ^ y.o.) 

U.h f. 

JiJ: j» 


where C 


’^'^ 1 * ••^ 8 '^ 1 " *^8 


is a numerical coefficient. 


The summation over integers i and ] is such that: 

i • 

I l'«l + I [jJ = > or 3 

«T| 

We will designate that: 

(r.6|= ^ ^ ... ^ 

and hence equation (7) takes on the form: 


(7) 


V.» 


We also make: 


10 
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£= + I 


£ = — I 


if 

if 


I C 


■ 



I 


I »«+ y ',= - • - 

M* 1 • 


Equation 


(6) is then written: 


dK 

dt 


i*m %.0 


( 8 ) 


Substituting the Lagrangian 
dk^/dt similarly yields; 

dk 


solution into the equation in 


Yf- = - I 


fc,.} - y £ X 




(?) 


We also calculate: 


dp, 
dt 

dq, 

dt 


- S [“•'•]{ 2/1,. q,-2/l„.q,.}+ y 


*.» 


y [•'• ' J|2/l,.;j,- 2.4„p,} - X »• X ^ „sin(v.th. 

v.» 


Short Period Terms 


Substituting the Lagrangian solution into the short period part 
of the perturbation function yields only short period terms that are 
first order with respect to the masses. It is only with the second 
mass order that we come across long period terms again. 

This time we have to consider for each planet the complete 
system of Lagrange equations (4), which we will write for a planet of 
subscript u in the form; 
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) 



( 10 ) 


We determine the short period effects argument by argument. For 

a short period argument i>^ + j being given integers, the 

functions F have the form: 

f, cos(» + y'A, ) + f . sin (I -r.;/., > 

where F and F are polynomials in h , k , q , h, . k , p , q , whose 
c s uu uvv v v 

coefficients are functions of (Iri the special 

case in which one of the two integers i, j is zero, i.e., in the case 
in which the short period argument takes on the form the 

functions F and F depend on h , k , p , q , n for v = 1, 2, 8 

and on the seven quantities for u ^ v.) 

We therefore substitute the Lagrangian solution into equations 
(10) , which after integration yield a short period increase in the 
Lagrangian solution. Then by doing a Taylor expansion of the second 
members of equations (10) , we obtain second order terms vrith respect 
to mass after substituting the first order that we have just found. 

We will retain only the second terms' long period parts. 

Since the masses are always substituted for numerically, the 
terms thus found in the second members of the Lagrangian equations 
have the same form as those coming directly from the perturbation 
function's fourth order long period terms. 

In contrast to the case of the perturbation function's long 
periods, for which we kept the fourth period terms, the criterion 
for choosing short period terms is numerical. What we did was to 
retain the beginning of the h, k, p, q expansion of all the arguments 
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causing changes in the Lagrangian solution frequencies of more than 
lO"^ "/yr. 


RESOLUTION OF THE SYSTE.MS OF DIFFERENTIAL EQUATIONS 


We saw that the contribution of the short period terms took on 
the same form as the terms coming directly from the perturbation 
function. We therefore have to resolve a system of differential 
equations having the form: 


dt 

dt 



(8) 

r*n W-8 

- z «-’] ( 2 > 1 .. /». + B,. /i, } -Y^^ 

(9) 


r*ii t.F 


as well as a similar system for variables p and q . 


For that, we are going to use the Krylov-Bogolyubov method. This 
method consists of finding a solution of the form: 


with 


J* 1 v.» 

^■= Y l«yAfjCOSVy-»- A/j, ^ ^COS(V.f>l 


I 

I 




• CT» 

? 8 

Z /i,>NyCosfl^+ Z 


^ n 


(It 


(10 


Jr 


^ =s +C 

-r V. ^ . 


( 11 ) 


/148 


( 12 ) 


By differentiating system (11) and taking account of (12) , 
we obtain: 




13 


ORIGINAL PAGE IS 
OF POOR QUALFTY 


dji 

\u 

dk 


i'- ‘ * .» 

8 

~ " Z sin X + By) — T («/. s)\fl ^ , sind/, B) 


y* I 


where 




(i;. s) - + • •■ -♦- igg, +y,.^j +72^2 + ■•• + . 


Furthermore, equations (8) and (9) yield, by plugging in (11) 

t//? { ® 

-jf- = I X ;,yA/,COSV'y+ j; ^ „COS(V.0)| 

I ly*=l *,# J 

f 8 11 

+ fi.. i I COS Vy+ X •W:.*..COS(V\f;) + X ^Vv .COS(ty«.f)) 

^■y f f * 

- X -'^y, X ''•yy'/,sinv’y+ X 8sin(v’.fl) 

»*y I ly=l V.« 


dt 


+ B„ 


X ^yV/ySin^’y+ X 

y=i V.* 


I ^ 

f - 1 £ X ?,.^.»sm(v%fl). 

I i.» 


Hence, we have two expressions, (13) and (14), for dh /dk and 
dk^/dt. We make equal their parts that are of third order with 
respect to variables h, k, p, and q: 


X B j/.^j Mj cos y^Uhs)^L.,,.»cos{Y,0) 

y=i v.» 

= X |-''^y X ''F, ,,.,COS(l/-.B)+ B,, X «COS(V.B)l + X >^y «. 8COS(V'. 

1 I • V.« J V.» 

J 

X X (ff- 


8 

V 

J~l 


K.0 


~ X t“-‘] |- **yi X •'^y.v.»sin(v’,B) + B,, X 'W, ^ «sin(i^’. fl)l + V ex^, ,sin(v-B). 

I V.» v.» I v“ 


The method now consists of establishing argument by argument 
identities within each order. Two cases arise: 

1) The case in which the argument (., 9 ) is equal to ij; . , i.e. 
have : ^ 

8 8 

X l'-.l+ X L/ml= I • 

^ • 1 W * I 


( 13 ) 


(14) 


(15) 


we 
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Here, we once again come across the Lagrangian solution argu- 
ments, and establishing identities between coefficients makes it 
possible to determine the new frequency values. The solutions then 
are expressed in Fourier series of these new arguments. Establishing 
the identities yields: 


and 




= X [m. .Vf. . 


r 


For a given >i<j, subtraction of these two equations furnishes: 
'll iu,v]x2A„. + gMM:, -M^J^ V A/.,J = 0. 

\r*» ' ••• 


The fact that 


and (u,v]B are not zero means 
uv 


that M = M' , whatever u and j are. 

U,4^j U,l[j 

We can then write: 

[«^ i] X 2.4., - y,) - y a/,.,, = . 

In the first member of this expression, we once again come across 
matrix E of the Lagrangian system. Subtracting the eigenvalue g^ 
from the principal diagonal means that the M . values (u = 1, 2, 

U , j 

..., 8) will not be independent. We then let: 

M . =0 


This is an arbitrary step in the Krylov-Bogolyubov method. It 

reduces to changing variables over the integration constants M ^ , 

N.. The choice of M. =0 does not specify the solution but 
D D'-j 

imposes the choice of a certain type of expansion for the coefficients 
of arguments 


Having made this choice, we then have for each 0 ^ , j 

system of eight equations in eight unkowns: B. and K 

n u , 4^ j 


fixed, a 
{u ^ j) . 
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2) Case in which the argument (^,0) is random, i.e. such that: 


I IU+ i l/J = 3. 

■•I 


Establishing argument by argument identities in equations (15) 
yields : 


and 




By subtraction, we obtain M = eM' ^ whatever u and argu- 

ment (ii',0) is. 


We can then write: 


This time there is no arbitrary step and the resolution of the 
system of eight equations in eight unknowns gives for each argument 
the eight values (u = 1, 2, 8). 


j“ I 




We therefore have expansions of and : 

■ 

K= Y. >i,;A/ySinv’, + X ^ <sinlt,-.0l 

v.« 

8 

E ^.Af^cosv'^-»- E £AL ,.,coi(v.O) 

And similarly we find: 


E Y -V., ,sin(v.^>» 

>« 1 ^.0 

■ 

E t^.jN,cosej+ E'-^...«co>(vj>i. 


j* I 


v.« 


(16) 


(17) 
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RESULTS AND DETERMINATION OF THE CONSTANTS OF INTEGRATION 

In Tables 6 and 7 we give the integration constants and frequen- 
cies of the solutions obtained solely from the second order h, k, 
p, q long period terms. Comparison of tables 3 and 6 show how great 
the contribution of the perturbation function's fourth order long 
period terms is, especially for Mercury, Venus, Earth, and Mars. 

Table 7 contains the frequency modifications anc C^, as well as 

the frequencies' new values: g. =g. +B., s. =s. +C.. 

1 1 1 1 1 1 

Table 6 

Constants of Integration (According to the Solution Based 
on the Perturbation Function's Second and Fourth Order Long Periods) 


i 

H 

ft 

.V, 


, 

0.IS454!<67 

83 .32 24: 16 

0.05887664 

11 5l 44.62 

2 

O.Ol8fr427!j 

IVi 26 15.44 

0,00323843 

302 59 26.61 

3 

0.0I2O4I0I 

3 18 18 13.08 

0.0096 ■’327 

22x 33 14.44 

4 

0.063 M07J 

307 01 46.20 

0.03227762 

2*5 46 53.61 

5 

0.042974SS 

27 17 43.13 

0.01.3841.157 

If* O' 34.95 

6 

0.048427S2 

127 29 49.50 

0.00786457 

125 49 11.16 

7 

0.03210686 

lOO 44 43.97 

0.008801 19 

316 0016.25 

8 

0.00932.'<59 

64 54 28.20 

0.00592089 

201 1107.71 


[Commas in tabulated material are equivalent to decimal points.] 


Table 7 

Modifications and New Frequencies in "/yr 
(According to the Solution Based on the Perturbation Function's 
Second and Fourth Order Long Periods) 


1 

B. 

ft 

C. 


1 

-0.258373 

+ 5.202996 

-0.44.39sj 

- 5.643943 


-O.OOi3'2l 

+ 7..345860 

-0.22o5lU 

- 6.'91897 

3 

-0,130032 

+ 17.201263 

-0,152'ta 

- is. 898909 

4 

-0.16916S 

+ 17.835416 

-0.23i.362 

- l'.>6:4'3 

5 

+ 0.0 1 8087 

+ 3.729488 

0 

0 

6 

+ 0.322115 

+ 22.608667 

- 0.6O69O7 

- ;6..348083 


+ 0.078.361 

+ 2.780148 

-0.08.3;ii5 

- 2.987531 

8 

+ 0.009180 

* 0.642296 

-0.0091.36 

- 0.686656 


[Commas in tabulated material are equivalent to decimal points.] 
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Table 8 

Constants of Integration (Complete Solution) 


i 

Af. 

A 



1 

0.I779I6U 

87 03 02-.09 

0.059627'2 

l2 0> 17.|5 

2 

0.02KM749 

193 35 04.97 

0.0031.3.338 

.3t'5 l3 i7.|9 

3 

0.00988829 

3 19 43 16.64 

001002295 

24-< 0| .<9.07 

4 

0.061 15I7.< 

307 48 39.56 

0.03l302'9 

2" 5606.16 

5 

0.(U34l6l6 

28 .30 11.60 

0.01383939 

l(>6iWil.5" 

6 

0.04S 14727 

127 42 54.52 

0.00785328 

1 25 3.'< 34.23 

7 

0.031261.34 

1 14 46 3l.5.'> 

0.OO8S0OSS 

3l6 1' 55.94 

8 

0.008‘39|8| 

"2 05 25.0.3 

0.0058SK06 

20: r 15.59 


[Commas in tabulated material are equivalent to decimal points.] 

Table 9 

Modifications and New Frequencies in "/yr (Complete Solution) 


* 8. 9. C. i, 


1 

- 0.262 2‘«) 

T 5.1990*7 

0.4109*9 

- 5.6 a 93" 

2 

-0.000490 

* 7.346091 

-0.199640 

- f.*7l027 

3 

-0.110749 

+ 17.220546 

-0.083094 

- 15.529299 

4 

-0.147321 

+ 17.85726.3 

-0.18265S 

- 1 *.5 18769 

5 

+ 0.495804 

+ 4.207205 

0 

0 

6 

+ 3.930206 

+ 26.216758 

-0.525894 

- 26.267070 

7 

+ 0.363394 

+ 3.065181 

-0.095511 

- 2.999837 

8 

+ 0.034747 

+ 0.66786.3 

-0.0l.39ii 

- 0.691431 


[Commas in tabulated material are equivalent to decimal points.] 


Tables 8 and 9 give the constants of integration and the frequen- 
cies for the complete solutions, i.e. the solutions that take the 
short periods into consideration. By comparing tables 6 and 8, we 
can see that the integration constants are once more greatly altered. 
Comparison of Tables 7 and 9 show that while the short period terms 
hardly change the frequencies related to the inside planets, the g^ 
and g^ frequencies on the contrary are changed to a much greater 
extent by the short period terms than by the perturbation function's 
fourth order long period terms. 

The modification of the constants of integration originates in 
the magnitude of the nonlinear terms found in particular in the 


V 


expressions of the elements related to the inside planets. Of course, 
we began by calculating these terms by numerical substitution of the 
Lagrangian solution in Tables 4 and 5. We then determined an analytic 
form of the expressions found so as to calculate the new integration /153 
constants. With the help of this analytic form and by making a first 
order Taylor expansion about the first values of the integration con- 
stants, we obtained integration constants of sufficient accuracy after 
several iterations. 


g Table 10 

X ,M. X 10°. Lagrangian Solution Amplitudes 
il D 


i N. 

t 


3 

4 

5 

6 

7 

8 

McTCUIV 

I779l6t3 

-2554951 

145862 

- 142815 

2426473 

11310 

62320 

708 

Venus 

U 19702 

2104749 

- 1 193243 

1254088 

1632845 

- 55183 

61366 

1076 

Earth 

39S077 

1642622 

988829 

-1250360 

1632798 

246415 

64990 

1250 

MaR 

65007 

291694 

2823123 

6115173 

1879988 

1609341 

86222 

1998 

iupiler 

689 

1163 

89 

46 

4341616 

-1555604 

217360 

5809 

Saturn. 

615 

- 1200 

702 

703 

3421845 

4814727 

198300 

6557 

I'tanus 

266 

293 

41 

38 

-4407122 

- 181041 

3126134 

137298 

Sepiune 

3 

M 

3 

2 

161243 

- 13497 

- 337321 

899181 



Q 

Table 11 






“i: 

Nj X 10 . 

Lagrangian Solution Amplitudes 




1 

2 

3 

4 

5 

6 

7 

8 

Mercury 

5962772 

- 1109444 

167847 

■ 72261 

■ (383939 

■ 13894 

-166511 

- 72418 

Venus 

562458 

31555S 

-1100237 

- 426565 

1 383939 

6010 

- 95861 

- 66262 

Earth 

405197 

254218 

1C02295 

280886 

1383939 

140436 

- 86594 

- 64951 

Mars 

86031 

56602 

- 1471377 

3130279 

1383939 

482015 

- 62835 

- 61532 

Jupiter 

986 

408 

7 

109 

1383939 

-315287 

- 47866 

- 58541 

Saturn 

- 1262 

5'6 

197 

- 1138 

1383939 

785328 

- 39025 

- 56429 

Uranus 

1057 

297 

16 

107 

1383939 

- 34725 

880088 

54753 

*‘fptune 

27 

16 

1 

12 

1 383939 

- 3851 

- 103543 

588806 

Lastly, we 

give the 

totality 

of our 

solution 

in Tables 8 to 

13. 


/154 


Hence, Table 8 contains the 32 integration constants. Table 9 gives 


the B. and frequency modifications as well as the frequencies' new 
values: + C^. The amplitudes of the Lagrangian 

solution corresponding to the new constants are given in Table 10 for 


I 




\ . 
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O g 

X. .M. X 10 and in Table 11 for p . .N . x 10 . Finally, Tables 12 and 

13 contain the amplitudes M n „ of the higher order argu- 

ments (^, 0 ), When computing these terms, we retained only the argu- 

-4 

ments whose amplitudes are higher than 10 for the planets Mercury, 
Venus, Earth, and Mars, and lO'^ for Jupiter, Saturn, Uranus, and 
Neptune. The zeros found in Tables 12 and 13 are amplitudes less 
than the retained percisions. 


Table 12 


Vcunwnt 

Mercury 

N'enus 

Earth 

Mars 

Argument 

Mercury 

Venus 

Earth 

Man 

Vt 

iSl 

0 

0 

0 

V* - 9, - 9, 

0 

0 

0 

- 177 


0 

0 

0 

- 199 

V’« - 9, + 9^ 

0 

0 

0 

- 240 


0 

0 

0 

453 

V '4 “9, +9, 

0 

- 536 

276 

6177 

U-Ot +»♦ 

0 

0 

0 

287 

V'4-29. 

0 

0 

0 

176 

v» 

0 

981 

- 1064 

- 13379 

V’4 -Ot+ftt 

0 

0 

0 

370 

Vt + 0, —9, 

0 

0 

' 0 

762 

Vt 

- 120 

456 

- 503 

0 

rs-vj-v. 

0 

0 

0 

134 

Vt 9* - i\ 

0 

0 

0 

- 354 

n-2$, 

254 

0 

0 

0 

Vt + 9, - 9, 

- 519 

4126 

- 3367 

- 14852 

Vj-9, -9, 

- 114 

0 

0 

0 

V'* + 9, -9» 

0 

0 

0 

2 ii 

V5-9, +0, 

102 

0 

0 

0 

V 4 + V 5 -V 7 

0 

0 

0 

118 

V,-9, +9j 

- 275 

0 

0 

0 

-V’4 - Vt 

0 

0 

0 

- 286 

rs~9, +9^ 

150 

0 

0 

0 

Vj - 2Vt 

0 

465 

- 378 

- '26 

V, 

14547 

2918 

:;so 

1092 

V) ■- V'4 - V« 

0 

0 

0 

- 133 

^5 + 9, -9, 

- 2172 

0 

0 

0 

Vi -- V’4 ■*- Vo 

0 

0 

0 

184 

rs+fl, -9, 

0 

- 216 

IK2 

503 

Vj -0| f 9^ 

396 

0 

0 

0 


0 

0 

0 

105 

V'j-9, -t 9j 

312 

0 

0 

0 

2r,-v7 

- 118 

0 

0 

0 

V’j - + 9, 

- 237 

0 

0 

0 

V'4-V>J + Vj 

0 

0 

0 

- ISO 

V'l “ 

- 275 

2457 

- 2354 

8S64 

■^<- 0 , +9, 

1644 

0 

0 

0 

Vi 9* + 9, 

0 

0 

0 

148 

+9, 

- 598 

0 

0 

0 

V'j 

0 

0 

0 

- 323 
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Table 12 


Ar^umeni 

Mercury 

Venus 

Earth 

Mars 

V, + (K - 

0 

0 

0 

- 162 

’<‘.1 ~ ^4 

115 

- 901 

8.57 

- 1422 

V) + - 0, 

0 

0 

0 

102 

•f'j + V’j ~ V’l 

0 

0 

0 

137 

Vj + Vi- V* 

0 

0 

0 

- 144 

-Vj - l»'i 

0 

- 691 

599 

388 

V:-2«, 

- 374 

0 

0 

0 

-Oj 

175 

0 

0 

0 

Vj ~f^i +0j 

- 137 

0 

0 

0 

V; ~ + ^2 

- 4711 

- 261 

- 174 

0 

V: 

2947 

0 

- 151 

0 

Vi 

0 

111 

0 

0 


- 875 

0 

0 

0 

V : * 

105 

0 

0 

0 

V: + H, -O7 

355 

0 

0 

0 

V: + Vj - Vi 

205 

0 

0 

0 

-2vj 

- 361 

0 

0 

0 

V1-V2-V5 

- 3340 

- I3l 

0 

0 

Vi -V: + Vs 

- 194 

0 

0 

0 

Vi - 'Cj + V* 

- 591 

0 

0 

0 

V.-2P, 

2231 

0 

0 

0 

Vi - 

- 1022 

0 

0 

0 

Vi — ~ 

350 

0 

0 

0 

V.-O, -«»i 

147 

0 

0 

0 

Vi -0| +0-, 

399 

0 

0 

0 

V, -«, +y. 

- 146 

0 

0 

0 

Vi-01 +0* 

2l2 

0 

0 

0 

V|-», +0J 

480 

0 

0 

0 

V, - fl, + O2 

- 8302 

0 

0 

0 

Vi-2«>2 

119 

0 

0 

0 

Vl - ®2 + 

- 139 

0 

0 

0 

Vi - flj +^i 

- 605 

0 

0 

0 

Vi 

0 

- 812 

- 650 

0 

V 1 + 

603 

0 

0 

0 

Vl +^2 -flj 

147 

0 

0 

190 

V, -»2 

67S6 

1372 

1007 

173 

Vl -fl, -fl, 

- 505 

164 

- 124 

135 

Vl + - ^*4 

- 209 

0 

0 

- 1087 

Vl -0, 

146 

0 

0 

0 

Vl +<*1 -^7 

- 316 

0 

0 

0 

Vl +Vj- V i 

645 

0 

0 

0 

Vl -V2- Vs 

258 

0 

0 

0 

2Vi - V: 

1544 

0 

.4^- 

0 

2Vi - Vs 

- 109 

0 

0 

0 

2Vi - 'Vi 

107 

0 

0 

0 

^V, - Vs 

- 824 

- 438 

- 337 

0 • 

Argument 

Jupiter 

Saturn 

Lranus 

Nep- 

tune 

... rt. -0, 

y 

_ 


5 

Vs 

0 

3 

2-1 

0 

+ rt- - 0, 

0 

0 

1 

1 

V- - 2V| 

0 

0 

0 

1 

V- -W 5 +^7 

0 

0 

t 

0 

V- - 0. 

0 

_ y 


0 

V- - 2'), 

0 

0 

4 

1 

- 0, 

0 

0 

4 

y 

..■■ - II. i (), 

0 

0 

S 

4 

V- - 2II, 

0 

0 

1 

0 


(cont. ) 


Argument 

Jupiter 

Saturn 

Uranuv 

Nep.. 

tune<( 

V- 

- 174 

- 198 

0 

- 

V- + (»- -f). 

0 

0 

5 

92 

V- + (V - »» 

0 

0 

- 2 

0 

V- +('. - "■ 

0 

0 

- 3 

0 

2v- Vs 

0 

0 

- 4 

3 

Vc,-2v- 

0 

0 

- 7 

0 

VV * V‘ * '/'a 

0 

1 

0 

0 

V., - - 1>„ 

0 

5 

0 

c 

V, - 2'V 

9 

15 

0 

0 

vv-'V - 

0 

0 

1 

0 

Vft - f*. - 07 

0 

0 

1 

0 

VV 20. 

0 

0 

1 

c 

VV-0- -0. 

0 

1 

0 

0 

v» 

- 594 

0 

192 

33 

V^^«^ -0- 

1 

0 

22 

3 

V6 - 0. 

0 

0 

1 

0 

v „ -r 0J - 0^ 

0 

- 5 

0 

0 

VV + V--V. 

0 

1 

0 

0 

2v,. - V* 

23 

78 

y 

0 

2vv - V, 

0 

3 

0 

0 

Vs-2v^ 

298 

- 962 

27 

0 

Vs- Vt-V- 

- 113 

348 

- 49 

0 

Vs - - Vl 

1 

4 

1 

0 

Vs - - Vs 

0 

0 

1 

0 

Vs - Vfr - V- 

0 

- 6 

- 36 

0 

Vs - 2v- 

7 

7 

- 440 

34 

Vs-V- - Vs 

0 

0 

- 3 

- 3.' 

Vs — V- ■ Vs 

1 

1 

- 13 

- 33 

Vs - 2Vs 

0 

0 

1 

1 j 

Vs -0s -0- 

0 

0 

- 5 

0 

Vs - 0* + 06 

1 

3 

0 

0 I 

Vs - 20, 

7 

5 

0 

0 . 

Vs ■ 06 - 0- 

0 

0 

2 

0 ' 

Vs - 0» - 0- 

1 

6 

- 5 

0 : 

Vs -20- 

0 

1 

- 8 

• 

Vs-0- -0s 

0 

0 

4 

2 

Vs -0- -0, 

0 

0 

3 

4 

Vs - 20s 

0 

0 

_ 2 

0 

Vs 

0 

257 

5257 

463 

Vs ■'■0- - 01 

1 

1 

15 

^1 

V, + 06 - 0- 

0 

0 

5 

0 

If', + 0, - 0, 

1 

- 3 

0 

0 

Vs 0s - 0- 

0 

0 

5 

0 

Vs T 0j - 0^ 

0 

0 

1 

0 

Vs + V- - V, 

0 

0 

2 

- i 

Vs ^ V, - V- 

III 

349 

15 
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CONCLUSION 


In this study of the long period variations of the olanetary 
elements, we added to the Lagrangian solution the terms of third 
order excentricity and inclination arising from the long period 
portion of the perturbation function calculated for the planets as 
a whole. We also took into consideration the influence of short 
period terms of second order mass. We particularly concentrated on 
determining the integration c;_>nstants that make the solutions agree 
with the mean elements when t = 1850.0 is used as time zero. 

The terms calculated with these constants are grouped together 
in Tables 8 to 13. Notice in these results the very strong coupling 
that exists, for a long period problem, in the planetary system. The 
magnitude of the terms arising from the short periods shov/s that 
there is no point to extending a theory to the fifth order on the 
basis of the perturbation function's long periods if the short periods 
are not taken into account. 

This work's essential task was therefore the comparison of the 
various effects according tf their origin so as to have an overall 
view of this problem and to be able to embark on the complete construc- 
tion of a long period theory. Our solution is in fact still incom- 
plete. Even so, we should take into account the direct terms of fifth 
order that must have an influence, especially forMercury, Venus, 

Earth, and Mars. We have yet to calculate the influence of short 
periods of higher orders of excentricity and inclination, and maybe 
even part of the third order with respect to masses in the case of 
the resonant argument 2X . - S'- between Juniter and Saturn. Such an 
investigation would be very important. Hov;ever, since the largest 
contributions have already been considered, it would no longer present 
any great difficulties for the determination of integration constants. 
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